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WALL BOUNDARY CONDITIONS. DISORDERED PHASE 
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Abstract. The six-vertex model, or the square ice model, with domain wall boundary 
conditions (DWBC) has been introduced and solved for finite N by Korepin and Izergin. 
The solution is based on the Yang-Baxter equations and it represents the free energy in 
terms of an iV x Hankel determinant. Paul Zinn- Justin observed that the Izergin-Korepin 
formula can be re-expressed in terms of the partition function of a random matrix model 
with a nonpolynomial interaction. We use this observation to obtain the large A'^ asymptotics 
of the six-vertex model with DWBC in the disordered phase. The solution is based on the 
Riemann-Hilbert approach and the Deift-Zhou nonlinear steepest descent method. As was 
noticed by Kuperberg, the problem of enumeration of alternating sign matrices (the ASM 
problem) is a special case of the the six-vertex model. We compare the obtained exact 
solution of the six- vertex model with known exact results for the 1,2, and 3 enumerations 
of ASMs, and also with the exact solution on the so-called free fermion line. We prove the 
conjecture of Zinn- Justin that the partition function of the six-vertex model with DWBC has 
the asymptotics, Zj\i ~ CN'^e^ ^ as A" ^ oo, and we find the exact value of the exponent 



1. Introduction 



The six-vertex model, or the model of two-dimensional ice, is stated on a square lattice 
with arrows on edges. The arrows obey the rule that at every vertex there are two arrows 

pointing in and two arrows pointing out. Such rule is sometimes called the ice-rule. There 
are only six possible configurations of arrows at each vertex, hence the name of the model, 
see Fig. 1. 

We will consider the domain wall boundary conditions (DWBC), in which the arrows on 
the upper and lower boundaries point in the square, and the ones on the left and right 
boundaries point out. One possible configuration with DWBC on the 4x4 lattice is shown 
on Fig. 2. 

The name of the square ice comes from the two-dimensional arrangement of water molecules, 
if20, with oxygen atoms at the vertices of the lattice and one hydrogen atom between each 
pair of adjacent oxygen atoms. We place an arrow in the direction from a hydrogen atom 
toward an oxygen atom if there is a bond between them. Thus, as we already noticed before, 
there are two inbound and two outbound arrows at each vertex. 
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(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



Figure 1 . The six arrow configurations allowed at a vertex. 



Figure 2. An example of 4 x 4 configuration. 



For each possible vertex state we assign a weight Wi, i = 1, . . . , 6, and define, as usual, 
the partition function, as a sum over all possible arrow configurations of the product of the 
vertex weights, 

Zn^ E ^<^ (1-1) 

arrow configurations 

where Ui is the number of vertices in the state i in a given arrow configuration. We will 
consider the case, when the weights are invariant under the simultaneous reversal of all 

a := W3 = W4, b := W5 = Wq, c := Wi = W2. (1.2) 
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Figure 3. The corresponding ice model. 



Define the parameter A as 

There are three physical phases for the six- vertex model: the ferroelectric phase, A > 1; 
the anti-ferroelectric phase, A < —1; and, the disordered phase, — 1 < A < 1. The phase 
diagram of the model is shown on Fig. 4. 




1 a/c 



Figure 4. The phase diagram of the model, where F, AF and D mark ferro- 
electric, antifcrroelectric, and disordered phases, respectively. The circular arc 
corresponds to the so-called "free fermion" line, where A = 0, and the three 
dots correspond to 1-, 2-, and 3-enumeration of alternating sign matrices. 
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In these phases we parametrize the weights in the standard way: for the ferroelectric 
phase, 

a = sinh(t — 7), 6 = sinh(t + 7), c = sinh(27), I7I < t, (1.4) 
for the anti-ferroelectric phase, 

a = sinh(7 — if:), 6 = sinh(7 + t), c = sinh(27), |t| < 7, (1.5) 
and for the disordered phase, 

a = sin(7 — t), 6 = sin(7 + t), c = sin(27), |t| < 7. (1.6) 

Here we will discuss the disordered phase, and we will use parametrization (|1.6|) . 

A solution for the free energy of the six-vertex model with periodic boundary conditions 
(PBC) was found by Lieb PZI^EHI by means of the Bethe Ansatz. In the most general form 
of the six-vertex model the Bethe Ansatz solution with PBC was obtained by Sutherland 
A detailed classification of the phases of the model is given in the review paper of 
Lieb and Wu f2P; see also the book of Baxter [ij. The six- vertex model with antiperiodic 
boundary conditions is solved in 

The six- vertex model with DWBC was introduced by Korepin in who derived an 
important recursion relation for the partition function of the model. This lead to a beautiful 
determinantal formula of Izergin fUj^, for the partition function of the six- vertex model with 
DWBC. A detailed proof of this formula and its generalizations are given in the paper of 
Izergin, Coker, and Korepin jT2]- When the weights are parameterized according to ()1.6|) . 
the formula of Izergin is 

[sin(7 + t) sin(7 - t)]^' 



T„, (1,7) 



where tat is the Hankel determinant. 



d 



■i+k-2. 



and 



\ / l<i,k<N 

m = , , tT^\ . (1.9) 

sm(7 + t) sm(7 — t) 

An elegant derivation of the Izergin determinantal formula from the Yang-Baxter equations 
is given in the papers of Korepin and Zinn- Justin T5j and Kuperberg [16j. 

One of the applications of the determinantal formula is that it implies that the partition 
function tn solves the Toda equation, 

tnTn -t'n = Tn+iTn-I, N >l, (') = —, (1.10) 

cf. 121] ■ This was used by Korepin and Zinn- Justin ^3] to derive the free energy of the 
six-vertex model with DWBC, assuming some Ansatz on the behavior of subdominant terms 
in the large asymptotics of the free energy. 

Another application of the Izergin determinantal formula is that rjy can be expressed in 
terms of a partition function of a random matrix model. The relation to the random matrix 
model was obtained and used by Zinn- Justin [2H]- This relation will be very important for us. 
It can be derived as follows. For the evaluation of the Hankel determinant, it is convenient 
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to use the integral representation of namely, to write it in the form of the Laplace 
transform, 



oo 

t\ 



where 



Then 



0(t) = / e*^m(A)rfA, (1.11) 



sinh§(7r-27) 

MA) = (1.12) 

smh 2 TT 



^ = / Ye*^m(A)dA, (1.13) 



and by substituting this into the Hankel determinant, ()1.8j) . we obtain that 

N 



TN= [ n[e*^'m(A,)rfA,]det(A^+'=-2)i<,,, 
^ i=i 

. N N 



;i.l4) 



i=l i=l 

Consider any permutation a G of variables Aj. From the last equation we have that 

^ N N 

TN= n[^*''^^(^^)^^^K-l)"det(A,'=-i)i<,,fc<ArJ]A^7,^. (1.15) 

i=l i=l 

By summing over a G Sn, we obtain that 

1 r ^ 
■ i=i 

(see j2HI), where A(A) is the Vandermonde determinant, 

A(A) = det(Af-^)i<,,fc<^ = HiXk - A,). (1.17) 



i<k 



mix] 



-Vix) 



Equation p.l6|) expresses rjv in terms of a matrix model integral. Namely, if 
then 

= jk^, J dMe^^^^''-y^^^\ (1.18) 

where the integration is over the space of N x N Hermitian matrices. The matrix model 
integral can be solved, furthermore, in terms of orthogonal polynomials. 

Introduce monic polynomials Pn{x) = + . . . orthogonal on the line with respect to the 
weight e*^m(x), so that 

/oo 
Pn{x)Pra{x)e*''m{x)dx = hJnm- (1-19) 
-oo 

Then it follows from (jl.l6j) that 

N-l 

= n (1-20) 

n=0 
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The orthogonal polynomials satisfy the three term recurrent relation, 

xPn{x) = Pn+l{x) + QnPn{x) + RnPn-l{x), 

where Rn can be found as 



Rn 



hn~l 



see, e.g., This gives that 



hn — ho Y\ -Rj) 
i=i 



where 



hn 



e ^m[x)dx 



sin(27) 



sin(7 + t) sin(7 — t) 
By substituting ()1.2H|1 into ()1.2()|1 . we obtain that 



N-l 



c n 



ra=l 



We will prove the following asymptotics of the recurrent coefficients Rr, 
Theorem 1.1. As n ^ oo, 

i2 



Rn 



n 



r 



[n 



-2-£\ 



R + cos{nuj) ^2 9"- + cn ^ + 0{ 

j: Kj<2 

where the sum is finite and it goes over j = 1, 2, . . . such that kj < 2, 

2 



R 



TT 



2cosf 



c 



t 

7 



e>0, 



2j 



uj = 7v{l + C); S = 1 + ^ _ 1 > 

27 



and 



where 



and 



where 



Also, 



276' 



COS 



•1) sm; 2^ , 



TTJ 



— - 1 

27 



TT 



2 



2^ 





2 








[/ 




TT 


yo 


vr 




vr 


2^" 


-( 


2^ 




7r7^ 



arg(/i + iy)f{fi)dfi + y\ny-y 



vr 



27 



sgn /i. 



vr" 



6(vr - 27) cos2 ^ 48cos2f 



;i.21) 
;i.22) 

1.23) 

;i.24) 

;i.25) 



;i.26) 



(1.27) 

(1.28) 
(1.29) 

(1.30) 
(1.31) 
(1.32) 
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The error term in ^1.2f^} is uniform on any compact subset of the set 

(t,7): |t|<7, 0<7<^} . (1.33) 

Remark. The method of the proof allows an extension of formula ()1.26p to an asymptotic 
series in negative powers of n. We stopped at terms of the order of n"^, because for higher 
order terms formula for Cj becomes complex. 

Denote 

^"^^'"(HS?- 

From Theorem II. II we will derive the following result. 
Theorem 1.2. .4s ^ oo, 

FN = F + OiN-^), (1.35) 

where 

F=-ln4 = ln (1.36) 

2 7' 27cosf ^ ^ 

This coincides with the formula of work 28j, obtained in the saddle-point approximation. 
Earlier it was derived in work [14j, from some Ansatz for the free energy asymptotics. For 
the partition function Z^v in (jl.7p we obtain from Theorem 11.21 the formula, 

i,.,>Z„^/ + 0(7V-') ,^,,( Mcos(20-c.s(2.)| \ 

\ 47 COS J 

Let us compare this formula and asymptotics p.26|) with known exact results. 
The free fermion line, 7 = f, \t\ < j. In this case the exact result is 

= 1, (1.38) 

see, e.g., [7j, which implies / = 0. This agrees with formula (jl.37p . which also gives / = 
when 7 = f . Moreover, the orthogonal polynomials in this case are the Meixner-Pollaczek 
polynomials, for which 

fl„ = J^ = !!!^. (1,39) 

cos"^ 2t 7^ 

cf. j7j. Thus, formula p.26|) is exact on the free fermion line, with no error term. This 
agrees with Theorem because from p.28j) , ()1.32|1 , cj = c = when 7 = f . 
The ASM (ice) point, 7 = | , t = 0. In this case we obtain from ()1.6|) that 

a = = c = — 
2 

hence 



;i.4o) 




;i.4i) 



where A{N) is the number of configurations in the six-vertex model with DWBC There is 
a one-to-one correspondence between the set of configurations in the six-vertex model with 
DWBC and the set oi N x N alternating sign matrices. By definition, an alternating sign 
matrix (ASM) is a matrix with the following properties: 
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• all entries of the matrix are —1, 0, 1; 

• if we look at the sequence of (— l)'s and I's, they are alternating along any row and 
any column; 

• the sum of entries is equal to 1 along any row and any column. 

The above correspondence is established as follows: given a configuration of arrows on edges, 
we assign (—1) to any vertex of type (1) on Fig. 1, 1 to any vertex of type (2), and to 
any vertex of other types. Then the configuration on the vertices gives rise to an ASM, and 
this correspondence is one-to-one. For instance. Fig. 5 shows the ASM corresponding to the 
configuration of arrows on Fig. 2. 



10 0^ 

1 

1-10 1 

10 0) 



Figure 5. ASM for the configuration of Fig. 2. 



For the number of ASMs there is an exact formula: 

This formula was conjectured in |22]; 12^]; and proved by Zeilberger |27|j by combinatorial 
arguments. Another proof was given by Kuperberg ^B]; who used formula (|1.7|) . The 
relation to classical orthogonal polynomials was found by Colomo and Pronko [7j, who used 
this relation to give a new proof of the ASM conjecture. The orthogonal polynomials in this 
case are the continuous Hahn polynomials and from we find that 

i^n = 4;^ = - + T^ + 0(n-)- (1-43) 



Formula ()1.2fi|l gives 



7r2 



' Oin 2 ^] 



4 144^2 

which agrees with p.43p . From p.42|) we find, see Appendix ^ that as — ^ oo. 



1.44) 



A(N)^C\m 7V-A(i-^ + 0(7V-»)). (1,45) 



where C > is a constant, so that 



^N = C[^X N-^^i^-J^^.+0{N-')], N^oo. (1.46) 
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Formula ()1.37|) gives / = In | , which agrees with the last formula. 
The X = 3 ASM point, 7 = | , t = 0. Here the exact result is 

3iV/2 

-A{N;3), 



;i.47) 



where 



A{2m + 1; 3) = 3'"('"+^) JJ 



k=l 



(3fc-l)! 
{m + k)\ 



A{2m + 2; 3) = 3m (3m + 2)!m' ^ ^_ _ 



:i.48) 



[(2m+ 1)!]2- 

In this case A{N\ 3) counts the number of alternating sign matrices with weight 3^^, where k 
is the number of (—1) entries. Formula ()1.48p for A[N\ 3) was conjectured in j^^l; and 
proved in JH] • The relation to classical orthogonal polynomials was again found by Colomo 
and Pronko [7^, who used it to give a new proof of formula ()1.48|) for the 3-enumeration of 
ASMs. The orthogonal polynomials in this case are expressed in terms of the continuous 
dual Hahn polynomials and from [7j we find that 



R 



'2m 



36m 



i?2m+i =4(3m + l)(3m + 2). (1.49) 

In this case the subdominant term in the asymptotics of i?„ exhibits a period 2 oscillation. 
Namely, we obtain from the last formula that 

-1 + i-lY 



Rn = + 



;i.5o) 



This perfectly fits to the frequency value = vr for C = in (|1.27|) . Moreover, formula (|1.26|) 
gives 



Rn 



36n^ 



+ 



-l)"Ci 



TC 



72n2 



+ 0{n 



1.51) 



which agrees with (jl.5U|) and it provides with the value of Ci = 
From ()1.48|1 we find, see Appendix El that as m — >■ 00, 



Ai2m;3)=Cs\^- 
where C3 > is a constant, and 



4m 



3-^(2771)^8 1 + 



77 



+ 0(A-3) 



A{2m+l-3) = C^[- 



so that 



and 



{2m+lf 



Ar2 



2m+l , X 1 / 131 



+ 0(777 



-3^ 



A{N;3) = Cs{-] 3--AT8 (1 + 



104 -27(-l)^ 
1944A2 



+ 0{N' 



' 1944A2 ^ 



N 



00. 



1.52) 



1.53) 



1.54) 



1.55) 



Formula (jl.37p gives / = In | , which agrees with the last formula. 
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We have the identity, 

Q'^Fn Rn 



(1.56) 



see, e.g., jS], which is equivalent to the Toda equation (jl.lOp . By using identity ()1.56|) . we 
obtain from Theorem 11.11 the following asymptotics. 

Theorem 1.3. As N ^ oo, 

^^:^^^Il = \cos{Nuj) J2 c.AT-. +ciV-2 + 0(iV-2-^). (1.57) 

This gives a quasiperiodic asymptotics, as — oo, of the second derivative of the sub- 
dominant terms. 

Zinn-Justin's conjecture. Paul Zinn- Justin conjectured in [28J that 

Z^r^CN^e^'^, (1.58) 

i.e., 

lim =1. (1.59) 

Formulae ()1.38p . ()1.46|) . and ()1.55|) confirm this conjecture, with the value of k, given as 

0, 7 = J, 1*1 <j; 

1 TT 

— , 7 = -, t = 0. 

I 18 ' 6 

Bogoliubov, Kitaev and Zvonarev obtained in ^ the asymptotics of Zj^ on the line ^ + ^ = 1, 
separating the disordered and antiferroelectric phases. This corresponds to the value 7 = 0. 
They found that in this case formula ()1.58|) holds with = j2 ■ 
With the help of Theorem II. II we will prove the following result. 

Theorem 1.4. We have that 

= CN^e^'^ (1 + 0{N'')) , e>0, (1.61) 

where 

1 272 ^ ^ 

K = -, 1.62 

12 37r(7r - 27) ' ^ ^ 

and C > is a constant. 

This proves the conjecture of Zinn- Justin, and it gives the exact value of the exponent k. 
Let us remark, that the presence of the power-like factor N'^ in the asymptotics of Zn in 
p. 6111 is rather unusual from the point of view of random matrix models. As was proven 
rigorously by Ercolani and McLaughlin ^0] , in the one-matrix model with an independent of 
N analytic interaction V{M) = M'^ + tVi{M), where t > is small, N'"^ In |^ is expanded 

into an asymptotic series in powers of A^"^. 

The set-up for the remainder of the paper is the following: 
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In Section 121 we describe a rescaling of the weight, which was introduced by Zinn- 
Justin [2H1, and which is convenient in the subsequent calculations. The rescaled 
random matrix model is described by a potential Vn{x) such that as — > oo, it 
approaches a limiting potential V{x). 

In Sections we evaluate the equilibrium measures for the random matrix models, 
first for the limiting one, corresponding to V{x), and then for the random matrix 
model, which corresponds to Vn{x). 

In Sectioniniwe remind the Riemann-Hilbert problem for orthogonal polynomials, and 
in Sections [Tlini we carry out the large asymptotic analysis of the Riemann-Hilbert 
problem, via a sequence of transformations and the Deift-Zhou nonlinear steepest 
descent method. 

We use the results of this analysis in Sectional where we obtain the large asymp- 
totics of the recurrent coefficients and prove Theorem ll.il The central point in the 
derivation of the subdominant asymptotics of the recurrent coefficient is a deforma- 
tion of the lenses boundary, see Section 11.11 During this deformation, the lenses 
boundary crosses poles of the function e~^'^^^^\ and every time it crosses a pole, a 
new subdominant term arises in the asymptotics of the recurrent coefficient. Section 
im gives a proof to Theorems I1.2H1.4I 

Finally, there are several Appendices to the paper, where some auxiliary results are 
proved and some exact large A^ asymptotics are obtained. 



2. Rescaling of the Weight 



Following ^8j, let us substitute Aj = in ()1.16|) . This reduces to 



where 



and 



tn 



N 

n 

i=l 



7 



7 



-i<C<i. 



The polynomials 

.N/ V 7 y 

are monic polynomials orthogonal with respect to the weight e^^^m ( — ) , so that 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



°° f Nx\ 

V 7 / 



where 



'y \ 2n+l 



N 



(2.5) 



(2.6) 
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It follows from (fT^ that 

N-l 

= Y\_ ^Nn- (2.7) 

n=0 

The polynomials PNn{x) satisfy the three term recurrent relation, 

xPNnix) = PN,n+lix) + QNnPNn{x) + RnuPn ,n-l{x) , (2.8) 

where ^ 

Rnu = ^n, Qnu = J^Qn- (2-9) 

In what follows we will evaluate the asymptotics of Rnn and Qn^ as N —>■ oo. In particular, 
for i^ATTv we will obtain the formula 

Rj^N = R + cos{Nuj) CjN-''^ +0{N-'^). (2.10) 

j: Kj<2 

Then ()2.9|) will provide us with the needed asymptotics of i?„ as n — > oo. 
Define 

^sinhiV/i(^-l 



sinhA^/if 



27 



Then 



hence 



(2.11) 



sinhiV/.^ V 7 ' ^ ^ 



AT 



/■ J][e-^^-(A'^)rf/i,]A(/i)2, (2.13) 
^ i=i 

Observe that as — > cxd, 

VNifi)^Vifi) = -Cfi+H (2.14) 

We will evaluate the equilibrium measures, first for V and then for Vat. But before we discuss 
some general formulae for equilibrium measures. 

3. Equilibrium Measure 

In this section we remind some facts concerning equilibrium measures, see jS], [HI- Let 
V{x) be a real analytic function such that 

lini ^ = 00. (3.1) 

ir^itoo In |x| 

The equilibrium measure, u^'^ = u^, for V is defined as a minimizer of the functional 

Iy(i,) = - \n\x - y\di^{x)du{y) + V{x)du{x), (3.2) 



over all probability measures z/ on M^. The minimizer exists and it is unique. The equilibrium 
measure has the following properties: 

• It is absolutely continuous with respect to the Lebesgue measure, di'^'^{x) = p{x)dx. 

• It is supported by a finite number of disjoint intervals, S = Pi]- 
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On S, 

1 " 

p{x) = —h{x)^/R{x + tO), Riz) = Y\{z-ai){z~pi), (3.3) 



1 = 1 



where h{x) is a real analytic function on the real line, and a/ R{z) is taken on the 
principal sheet. 

The function h{x) is expressed by the contour integral, 

over any closed contour F, with S in its interior. For the equilibrium measure, consider its 
resolvent, 

u{z)= .eC\^, (3.5) 

Js z — n 

and the (^-function, 

9{z)= p{fJ')^og{z-fi)dfi, z e C\{-oo,pq], (3.6) 
Js 

where for log z the principal branch is taken. Then 

g'{z)=u;{z), (3.7) 
and, by the jump property of the Cauchy integral, 

+ iO) — Lj{fi — iO) = —2'Kip{fi), p E S. (3.8) 

As 2 — > oo, 

uj{z) = z-^ + Oiz-^), giz) = \ogz + Oiz-^). (3.9) 

The equilibrium measure is uniquely determined by the condition that there exists a real 
constant / such that 

• For any p E S, 

g{fi + iO) + g{p~iO) -V{n) -I = 0, fi e S. (3.10) 

• For any p E'R \ S, 

g{fi + iO) + g{fi~iO) -V{fi) -I <0, peR\S, (3.11) 
see jH], [SI- It implies the equation, 

uj{p + iO) +uj{fi-iO) = V'{p), fxES. (3.12) 
A solution to this equation can be found as 

^' 2m Js{z-x)^R{x + iO) 
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2'Ki I p{s)xs{s)ds, ai < /i < (3q, 



g{fi + iO) - g{n-iO) = { 2ni I p{s)xsis)ds, ai < fx < pg, (3.14) 

where xs is the characteristic function of S. 

We will be interested in the case when V is convex. In this case the equilibrium measure 
is supported by one interval, say, and ()3.14p reduces to 



2'Ki / p{s)ds, OL < p < (3, 
Jfi 

0, p>f3. 



(3.15) 



For zo E C and r > 0, we will us the standard notation for the disk, 

D{zq, r) = {z eC: \z-Zo\<r}. (3.16) 
From ()3.3|) . ()3.10|) and ()3.15|) . one obtains that there exists r > such that 

-2g{z) + V{z) + l= [ h{s)^y{s-a){s-(3)ds, z e D{P,r) \[a, p], (3.17) 

and 

roc 

-2g{z) + V{z) + 1 = -2Kisgn(lmz) + h{s)^J{a - s){(3 - s)ds, z E D{a,r)\ 

J z 



Finally, it follows from ()3.9|) that the function e^^^^"^ is analytic on C \ and 



z ^ oo. 



(3.18) 



(3.19) 



4. Equilibrium Measure for V 
In this section we consider the equilibrium measure for the potential 

r(/i) = -C/i+|/i|, |CI<i. 

It is obviously a convex function, and 

V\p) = -C + sgn(/i). 
Integral ()3.13p is explicitly evaluated in this 



1 — C 2 
uj{z) = — — + — log 

2 ITT 



^/ (3{z — a) — i^ —a{z — (3) 



a 



and from the asymptotics. 



uj(z) = z'^ + Oiz 



z ^ oo, 



one finds that 



TT 71 

a = — TT tan — (1 — C), = tt tan — (1 + C), 



(4.1) 
(4.2) 

(4.3) 

(4.4) 
(4.5) 
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see |2H1- Observe that 



j3 + a = 27rtan 



[3 — a 



277 



COS -Y 



{-a)(3 = ttI 



(4.6) 



For the square root in ()4.3p we take the principal branch, with a cut on the negative half-axis. 
The function uj{z) is analytic on C \ [a, (3]. On 



UJ[II±lO) = ± — l0£ 



a < fi < P, (4.7) 



so that the density function p{fi) is equal to 



Pif^) = log 



A//5(/i - a) + ^/-a{p - n) 



a < n < 13. 



Observe that p(/i) has a logarithmic singularity at the origin. From ()4.7j) . 

-C-1 -C + 1 



a;(a) 

From (p3|l we obtain that 



^(z) = zcj(2)+21og (yF^ + - /^) -(1 + 2 log 2) 
(see Appendix iBl below) . This implies that 



^ 2 
p{s)ds = —fip{fi) H — arctan - 



TT 



fj, — 



a<fi< p. 



In particular, 



^ 1 + C 

p{s)ds = 



In addition, we have that 

g{p + iO) + g{p-iO) -V{p) -I = 0, a < p < (3. 
By (gH), uj{P) = i^, hence 

I = 2g{(3) -V{(3) =21n(/?-a) -2-41n2. 
Define an analytic continuation of the potential V{p) = —C,p +1/^1 from R to C as 

[l-Qz, Re^>0, 
'-l-Qz, Rez<0 



V{z) 



(4.8) 

(4.9) 
(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 

(4.15) 



(the function V{z) is two- valued on Re 2; = 0). In what follows we will use the function h{z) 
defined by the formulae. 



h{z) 



7c^{z -a){z - (3) 



log 



^yp{z -a) - iy/{-a){z - (3) 



a) 



Rez>0. (4.16) 
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h{z) 



Ai 



log 



TT 



^{-a){f3-z)+t^(3{a-z) 
y/i-z){P-a) 



Rez < 0, 



(4.17) 



where the square root and logarithm are taken on the principal sheet. The function h{z) has 
the following properties: 

(i) h{z) is analytic in C \ iM and 

h{a) = > 0, h{P) = -rr^ > 0; 



h'{a) 



-a) {P — a) 

4(13 -3a) 
Sa^ (/? — a) 



2 ' 



l3i/3-a) 

4(a - 3(3) 
3/32 {(3 -a 



>2' 



(ii) by gM, 

(iii) by (Q, 



uj{z) 



V'{z) h{z)^{z-a){z-P) 



(4.18) 
(4.19) 

(4.20) 

(4.21) 



2 2 

Asymptotic formulae for orthogonal polynomials with weight 1)4.11) with C = were obtained, 
via the Riemann-Hilbert approach, in paper ^Hj by Kriecherbauer and McLaughlin. In fact, 
they studied a more general case, of the Freud potentials of the form V{fi) = 

5. Equilibrium Measure for Vn 
From ()2.11|) we obtain that 



(- 

V27 



1 cothA^/i 



(- 

V27 



— cothA^u C- 

27 ^^27 ^ 



(5.1) 



The function V^(/i) is increasing, hence V]\f{fi) is convex. Its equilibrium measure, pN{fJ')dfi, 
is supported by one interval [a^, (3n]- As A^ — 00, the equilibrium measure for Vn converges 
to the one for V. In this section we will derive some asymptotic formulas for a^, I3n and 
PAr(/i) as A^ ^ 00. Consider the resolvent. 



ujn{,z) 



z — n 



z eC \ [on^Pn]- 



Then 



and 



Pn{p) 



[u;Ar(/i + ZO) - to'Ar(/i - zO)], On < p < (3^ 



'N, 



f 

UJNifJ' ± ^0) = =F7r2p7v(/^) + P-V- 

J a 



PN{x)dx 



a]\f < p < P 



'N, 



(5.2) 

(5.3) 
(5.4) 



' ajv A'' 

where P.V. J is the principal value of the integral. The resolvent solves the equation, 

ujn{p + iO) + LUNip - iO) = Vl^{p), aN<p<PN- (5.5) 
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The solution to this equation is 



where 

i?iv(z) = (2; - aiv)(2; - /?7v), (5.7) 
and ^^Rn{z) is taken on the principal sheet, with a cut on [aNil^N]- As usual, 

\/ Rn(x), = lim \/ Rn(x + is). 

Evaluation of the end-points. From ()5.2j) we have that ujm{z) ~ ^ as 2; 00. By eval- 
uating the large z asymptotics of the integral on the right in ()5.6|) . we obtain the equations, 

^ ^^^''^ rfx = 0, (5.8) 



2vr A/(a; - aAr)(/5Ar - x) 
and 

J r ^ -'^^'-> - (5.9) 

2^ ia^^ V (a; - Otv) (/^AT - x) 

From these two equations we obtain the following asymptotics of oat, (3^ as — >■ 00. 
Proposition 5.1. As N 00, 



,72 (2 sin ^-1) ^ 

3(7r-27)cosf ^ 

^72(2sin^ + l) , 
3(7r — 27) cos 



(5.10) 



where a, (3 are given in 

Proof of Proposition 15.11 is given in Appendix O below. 

Evaluation of the density. Consider now the asymptotics of the density function Pn{x)- 
As ^ cxD, it approaches the density function p{x) given in ()4.8|) . The density p(x) has 
a logarithmic singularity at x = 0. For the singularity is smoothed out and we are 
interested in the large N asymptotics of pat near the origin. From ()5.3|) and ()5.6|) we obtain 
that 

. s x/tnUA f^" Vlr(x)dx 

Pn{p) = -^h^Py- / 7 , , aM<p<PN. (5.11) 

^'^ JaM (p-a;)Vr7v(x) 

where 

r^ix) = (x — aN){PN — x). (5.12) 
Observe that Pn{p) is analytic for < /i < /^a?. From ()5.ip we have that 

VUx)=sgnx-C + f{Nx), (5.13) 

where 

fix) = — cothx 1 cothx 1 — sgnx, (5.14) 

^ ^ 27 27 V27 ; V27 ; ^ ' 
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27r2 



P.V. 



■Pn (^sgnx — Qdx 



'On 



(/i - x)^/rN{x) 



aN < p < Pn, 

OTV < /i < /^AT. 



The function p%{p) is evaluated exphcitly as 



71' 



\/Pn{p- aN) + \/^aN{(3N - p) 



\/\p\{I3n - oin] 



aN < p < P_ 



'N- 



[cf. (gSl)]. Set 



Then 



2 

+ — log 

2 ITT 



\/z{(3n - (^n) 
[cf. (Q]. 

To describe the large asymptotics of p^ip), introduce the function, 

From expilicit formula ()5.14|) . we have the following properties of the function f{x) 

• /(x) satisfies the estimate, 

|/(x)|<Coe-^N, 

with some Co, C > 0; 

• f{x) is an odd function; 

• the function f{x) + sgnx is real analytic. 

These properties of f{x) imply the following properties of the function k{p): 

• the function k{p) is even and 

k{p) = —2 In + kQ{p), 

where ko{p) is real analytic; 

• as /i — *• cxD, 



where 



kip) = ^ + 0{p 
p^ 



C = xf{x)dx 



p ^ oo, 
2717^ 



(5.15) 
(5.16) 

(5.17) 

(5.18) 

(5.19) 
(5.20) 



(5.21) 



(5.22) 



(5.23) 



(5.24) 



3(7r - 27) ■ 

We use the properties of k{p) to prove the following asymptotics of the function p\[{p) 



(5.25) 
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Proposition 5.2. As N oo, uniformly in the interval aN ^ Pn, 



27r2 



k{Nfi) + 0{N~ 



(5.26) 



In addition, there exists a family of complex domains {fir, r > 0} such that [a^+r, f3N—r] C 
Qj. and Qr C Qr' whenever r > r' , and such that the function 



(5.27) 
(5.28) 



can be analytically continued to Qj. and as N oo, 

sup \er,{z)\=0{N-'). 

The proof of Proposition 15.21 is given below in Appendix |Dl 
Proposition 15.21 implies that 

PNip) = P%ip) - ^ k{Nn) + 0{N-^), fie K,/?iv], 

and this equation can be extended to the complex domain n^,, r > 0. This can be further 
specified as follows. Let r > be an arbitrary fixed number such that r < ^mm{— a, j3}. 
Then 



(5.29) 



For /i G [a AT + r, (3n - r], 



Pn{p) = pip) 



2n' 



k{Nfi) + 0{N- 



where p(/i) is given in ()4.8|) . 
• For /i G [ajv, Otv + r] U [Pn - r, Pn], 

pM=p%{fi) + 0{N-'). 

Observe that ()5.30|) implies that 



Pn{p) = ^IriN + a{p) ^ u ^ at.a , at-2 



^ - ko{Np) + 0{N-'), /i G K + r, pN - r] 



where 



lo£ 



^Jl3{fi - a) + A/-a(/3 - ji) 



VP - " 



a < ji < p. 



(5.30) 

(5.31) 
(5.32) 
(5.33) 



From ()5.29|) we obtain the following result: 



Proposition 5.3. We have that 



PN{p)dp 



1 + C 



+ 0(iV- 



Proof. By an explicit integration of ()5.18|) we have that 



/3jv 2 

p%{x)dx = —pp%{p) H — arctan - 



Pn - P 

p — ttiv 



OlN <P<P. 



N- 



In particular, 



/3jv 



p%{x)dx = 1, 



(5.34) 



(5.35) 



(5.36) 
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and 

pPn 1 I /■ 

p%{x)dx = ^ + OiN-'). (5.37) 

Since 

/•/3jv 

/ pN{x)dx = 1, (5.38) 

we obtain from ()5.29|1 that 

/oo 
k{x)dx = 0. (5.39) 
■oo 

Since k{x) is even, this imphes that 

POO 

k{x)dx = 0. (5.40) 



Therefore, 



1 + C 



PnWp = / p'Mdp - — / kiNp)dp + 0(iV-2) = ^ + 0(iV-2). (5.41) 
Jo ^'^ Jo ^ 

Proposition 15.31 is proved. 

Evaluation of the resolvent. The large asymptotics of the function ujn{z) can be 
obtained as follows. By ()5.15|) and ()5.27|) . 

PN{p) = P%{p)-7^k{Np)+eN{p), (5.42) 

hence 

u;n{z) = uj%{z) - ^ m{Nz) + (5.43) 

where 

miz) = r = r (5.44) 

Observe that ^n{z) is an analytic function in {Qr \ M), r > 0. Consider a complex domain 
Ur such that the closure of Ur belongs to fir and [a^ + r, — r] C Ur- Then from (j5.28j) 
and the analyticity of e^lz) in Qr we obtain that 

sup \^^{z)\=0{N~^). (5.45) 

zec/r-\K 



We have that 

m{z) = sgn{lmz)7Ti f IM^!^^ (5.46) 
y_oo z-p 



Indeed, if we introduce the Fourier transform. 



k{T) = ^ I e~'^^'k{p)dp, (5.47) 



27r 

then ()5.44j) implies that 

m±{T) = T'2ni9{±T)k{T), (5.48) 

where m±{p) = m{p±iO), and 6{t) = 1 for r > and 6{t) = for r < 0. Also, from ()5.21|) . 

fc(r) = -7risgn(r)/(r), (5.49) 
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r 

m(/x ± iO) = ±TTi i 
J —I 



(5.50) 
(5.51) 



hence 

m±(r) = -27r^^(±r)/(r). 
By taking the inverse Fourier transforms, we obtain that 

f\x)dx 
/i ± iO — X ' 

which imphes, by means of analytic continuation, (j5.46|) . 

By using the hsted above properties of /(x), we obtain from ()5.46|) the following properties 
of m(z), z e (C\M): 

• The symmetry conditions, 

m{—z) = —m{z), m{z) = m{z). (5.52) 

• The representation, 

m{z) = ±2711 log z + niQ^z), ±Im2;>0, (5.53) 

with log 2 on the principal sheet, where mQ{z) is analytic in the closed half-planes 
{±lm z > 0}, and 

(5.54) 
(5.55) 



As z — i> oo, 



mo{-z) = -mo{z), mo{z) = mo{z). 
Csgn (Im^;) 



m[z) 



+ 0{z- 



We summarize the properties of ujn{z) in the following proposition. 

Proposition 5.4. For any r > there exists an independent of N complex neighborhood Ur 
of the interval aN + ^ ^ 13 n ~ such that for z ^ Ur, equation ^5.4^ holds, in which 
uj%{z) is given by \5.2(^) . m{z), by (|5.^6| ), and ^n{z) satisfies estimate Ii5.4^ - In addition, 
for z G Ur, 



ilnN 1 
uJn[z) = T^— + b{z) - —mo{Nz) + 0{N-^), ±Im 2; > 0, 



TT 



27r2 



where 



biz) 



2 

^— + — loj 

2 tTT 



a/ (3{z — a) — i^ —ai^z — (3) 



VP - « 



with a cut on [a, [3], and 

mQ{z) = ±7ri 
with log 2; on the principal sheet. 



fip)dfx 
z — n 



2 \ogz 



±\mz > 0, 



(5.56) 
(5.57) 

(5.58) 



Observe that both b{z) and mo{Nz) have a jump across [a,/?], and 



b{z) = b{z), mo{z) = mo{z). 



By using ()4.fi|l . we find that 



b{+iO) 



1-C , 1 



+ — In 27r cos — . 



(5.59) 
(5.60) 
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Evaluation of the constant of integration. Let us evaluate In- By for any 



Take = |. By (I2.11|l . there exists c > such that 



Also, by (ITHIl . 



9N- \ f ) + ( f 



2 / pAr(x)ln 



/5 



By ()5.29j) and ()5.1()j) . we can reduce this to 



9n~ 1 + gN+ = 2 f p{x) In 



f3 



dx 



1 



I3n 



k{Nx) In 



"AT 



/3 



— x 



dx + 0{N~ 



From (fKT^ and (jO^ we obtain that 



/3jv 



A;(A^x) In 



/5 

2-^ 



c/x = 0{N~ 



hence 



Thus, 



p(x) In 



2"^ 



dx + 0{N~ 



)+0(iV-2) = / + 0(iV-2), 



where by (imil . Z = 2\n{(3 - a) -2-41n2. 



(5.61) 



(5.62) 



(5.63) 



(5.64) 



(5.65) 



(5.66) 



(5.67) 



6. RiEMANN-HlLBERT PROBLEM 

The Riemann-Hilbert (RH) problem for orthogonal polynomials with respect to the weight 
w{ji) is the following: 

(i) (analyticity) Y{z) = (Yij{z))^ ^is a. matrix valued analytic function on C\R which 
has limits on the real line, Y±{fi), so that for all A > 0, 

(6.1) 



(ii) (jump condition) 



lim max \Y(jj,±ie) — Y±(jj,)\=0. 

e^0+ -A<ij,<A 



(6.2) 
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(iii) (asymptotics at infinity) 

Y{z)=[l + 0{\zr)](^'^ \z\^oo. (6.3) 

Proposition 6.1. The RH problem (i)-(iii) has a unique solution given by 



Y{z) — ( 27ri7r„,i(2) r -tt j-^i(^]^7M)'^M ) (^'^^ 

where vr„(/i) = /!"■ + ... denotes the n-th monic orthogonal polynomial with respect to the 
measure w{fi)dfj, and hn = J^7Tn{lJ'yw{ii)dfi. Furthermore, there exist 2x2 matrices Yj, j = 
1,2, ... , so that for all m > 1, 

Y(z) ° ) = / + ^ + If + . . . + + O (|.|-»-) , \A - (6.5) 



and 



27rz 

hn = -27li(Yi)i2, hn-l = 



(n)21 

Rn = (Vl)2l(n)l2, (6.6) 

Qn = TTpr^ + (^l)ll) 
l^lj21 

where Qn, Rn o'^e the recurrence coefficients associated to the orthogonal polynomials, 

ZTlniz) = Tln+liz) + QnT^n^z) + i?„7r„_i(z). (6.7) 

RH problem (i)-(iii) and Proposition 1 liold for a general weight w(yu) (see [2] for conditions 
on w{fi)). In our case 

w{fi) = e-^^^(^), (6.8) 

and ()6.6p reads 

27ri 



^ATn = -27ri(Yi)i2, hN,n-l= . 

Rnu = (n)2i(>l)i2, (6.9) 

Qnu = TTTT^ + C^l)!!) 
I^lj21 

7. Transformations of the RH Problem 

We will follow jH] to find the asymptotics of the solution Y{z) to the Riemann-Hilbert 
problem (i)-(iii) in the case when n = N and N ^ oo. 
Transformation of the RH problem (|6.1|) - (|6.3|) . Set 

T(z) ^ e-^^'^^F(;^)e-^M^)-*)'^^ zeC\R; = (^J . (7.1) 

where 

9n{z) = 1 pN{fJ')^og{z - fi)dfj,. (7.2) 
Then T{z) solves the following RH problem: 
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(i) (analyticity) T{z) = (Tij{z))^ 2 is a matrix valued analytic function on C\M which 
has limits on the real line, T-t(/i), so that for all A > 0, 

lim max |T(/i ± ie) — T±(/i)| = 0. (7.3) 

£^0+ -A<ii<A 

(ii) (jump condition) 

TM=T-{fi)Mfi), (7.4) 

where 

/^-i^{9N+M-gN-M) I \ 

= { Q g7V(gjv+{/^)-9^-{/.)) J ' e [q;7v,/37v], (7.5) 

and 



JT(/i)=fQ ^ j, /iGM\K,/5^]. (7.6) 

(iii) (asymptotics at infinity) 

T{z) = I + {\z\~^) , |;z|-^oo. (7.7) 

The key point here is that the (21) element of the matrix Jrif^) on [aNyf^^] is equal to 1, 
due to equation ()3.1()|1 . For convenience, let us rewrite the recurrent coefficients Qnn, Rnn 
in the new terms: 

fiNN = -27rze^'^(Ti)i2, hN,N-i = — 77^ — 

(J 1)21 

Rnn = (Ti)2i(Ti)i2, (7.8) 
Qnn = -Ltt^ + (^i)ii- 
Jump matrix factorization. Denote for the sake of brevity 

G'7v(/i) = gN+ifJ') - gN-{fi)- (7.9) 
There is the following factorization of the jump matrix Jt on [aN,l3N]'- 

e-^«-(^) 1 1 0\/0 A/ 1 0\_ 

(,^GMM]-{eNG^(p) iM-1 ] I e~^^^(^) ^]=V-Vov+. [(.iU) 

Substituting this factorization into ()7.4j) for fi G [aN,l3N], we obtain that 

T+{fi) =T_{fi)v-{fi)vov+{fi), (7.11) 

or 

[T+{ix)v-\f^)] = [T4fi)v4fi)]vo, fi e [aN,PN]. (7.12) 

Lenses. By using the factorization of jump matrix ()7.1()|1 above, we can transform the 
RH problem for T in the following way. Consider the contours S]^ and E]^ on the complex 
plane from aj\f to Pn, as shown on Figure 6. 
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'N 



Figure 6. The lenses. 



The contours go closer and closer to the origin as ^ oo. Namely, by ()2.12|1 



/ ^ sinhA^^TT- 



sinhA^z(^ - 1 



(7.13) 



so that the function e^^^^^^ has poles on the imaginary axis. Consider the first pole in the 
upper half-plane, 

^1 = ^^- (7.14) 

The contour should be in the upper half-plane and it should cross the imaginary axis 
below zi, say, at \zi. We take = E^. We call the region between E^ (respectively, E^) 
and [aN^f^N] the upper (respectively, lower) lens. Let 

{T{z), outside of the lenses, 
T{z)[v+{z)\~^ , in the upper lens, 
T{z)v^{z)^ in the lower lens. 

Then S{z) solves the following RH problem: 

(i) (analyticity) S{z) is analytic on C \ (M U E). 

(ii) (jump condition) 

S+{z) = S^{z)Js{z), 



(7.15) 



(7.16) 



where 



Js(z) 



l gA^[9]V+(2)+3JV-(z)-VAr(z)-'jv] 
1 



zeM\[aN,P_ 



1 
-1 



z e [ajv,/3^ 



^ ° 1 ^ ^ E+ 



(7.17) 



1 



z G E^. 
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(iii) (asymptotics at infinity) 

S{z) = I + 0{\z\-^) , 1^1 ^oo. 

Evaluation of the functions e=^^'^^(^). By and ^J^ . 

rf^N 

G'Ar(;U) = 2TTi / p]sf{s)ds, aN < /i < Pn, 
J II 

hence, in particular, by ()5.34|) . 

rf^N 

Gn{0) = 2m / pNis)ds = 7ri(l + C) + OiN~^). 
Jo 

Consider first e"^'^^*-^-' on S^. From ()H.14p we have that 

gN-ifJ') = Vn{p) + In- 5'iv+(yu), aN < fi'< Pn- 



(7.18) 



(7.19) 



(7.20) 



(7.21) 



The RHS of this equation is extended to Im /i > and this gives us an analytic continuation 
of 5'Ar_(/i). By applying this continuation to ()7.9p . we obtain that for z G S^, 



hence 
By (EH, 

In particular. 
Therefore, 

so that 



Gn{z) = 2gN{z) - Vn{z) - In, 



^-NGn{z) ^ ^-2NgN{z)+NVN{z)+NlN _ 



sinhiVz^f 

27 ^-N(z 



sinh Nz{^-1 



,NVn{0) _ !L 



-NGNiz) 
-NGn{0) 



TT — 27 



TT 



vr — 27 



-2N[g,^{z)-gM{0)]+NVM(z) 



NGn{z) ^ g-2Af[3iv(2)-Sjv(0)] 



sinh Nz 



^ ^-Niz 



sinh Nz i^-l 



where 



By (iinnD, 



TT 



= !^J2e-^-'(^+«)(l + O(iV-i)). 



TT 



Observe that if Im z > then 



gN{z) - gN{0) = / UN{s)ds 



(7.22) 
(7.23) 

(7.24) 

(7.25) 
(7.26) 

(7.27) 

(7.28) 
(7.29) 
(7.30) 
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where the integration is taken over the interval [0, z\. From ()5.5(')|) we obtain that 



9n{z) - 5'7v(0) 



iz hiN 



TC 



mo{Ns)ds + 0{N-^\z\). (7.31) 



In particular, for z = iN ^y, where y > is bounded, we obtain that 

1 



gN{iN-'y) - gM{Q) = + ih{+iQ)yN-^ 



Mo{iy)N'^ + 0{N 



-2\ 



where 



Thus, (fri7\i gives that 



mo{s)ds. 



sm y 



27 



(l + 0(iV-i)), y>0, 



where 
and 

with 



U = -71(1+0, 



kN{y) 



TT — 27 



TT 



if{y) = -i2b{+t0)y + ^Moiiy) - iCv- 



By using the value of b{+iO) given in ()5.60|) . we obtain that 



21/ 



TT 



In ( 27r cos 



;Mo{iy). 



From ()7.33|1 and ()5.58|1 we obtain that 

\Mo{iy)=Q{iy)-Q{+iQ), 

where 



Q{z) 



TT 



Since /(/i) is odd, we have that 



Q{ty) - Q(+zO) 



TC 



\og{z — n)f{ji)djj, — 2z log z + 2z 



w:g{iy + fi)f{fi)dfi + y\ny-y 



+ iy. 



Thus, 



^(y) = -— In (271 cos ^ 



2 

TC 



aTg{iy + fi)f{fi)dfi + ylny-y 



(7.32) 
(7.33) 

(7.34) 

(7.35) 
(7.36) 

(7.37) 
(7.38) 

(7.39) 
(7.40) 

(7.41) 
(7.42) 



Consider now Imz < 0. Similar to ()7.22|1 we have that G'Ar(/i) is analytically continued to 
G]\f{z) with Imz < as 



Gn{z) = -2gN{z) + Vn{z) + I 



N, 



Im z < 0. 



From ira . 



9Niz) = gNiz). 



(7.43) 
(7.44) 
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Also, Vm(z) = Vat (-2) and /at G M, hence 



Gm{z) = -Gn{z). (7.45) 

From fj7.:^4)l we obtain now that 



^NG^i-^N-^y) ^ ^-NG^{iN-^y) ^ e-'^^^k^iy) , (1 + 0{N-^)), y > 0. (7.46) 



Model RH problem. Note that the jump matrix Js{z) converges, as — *• cxo, to the 
identity matrix, except on the interval [a, /?] where it is constant. This leads to the following 
model RH problem. 

(i) M{z) is analytic on C \ [a,/?]. 

(ii) (jump condition) 

M+{z) = M4z)Jm, z G K,/5^], (7.47) 

where 

Jm=(^, iV (7.48) 



(iii) (asymptotics at infinity) 

M{z) = I + {\z\~^) , l^l^oo. (7.49) 

Solution to the model RH problem. The model RH problem can be solved explicitly. 
Namely, let us reduce it to a pair of scalar RH problems that are solved by the Plemelj- 
Sohotski formula. By diagonalizing the matrix Jm-, we have that 

A _ 1 A 1 \ A Q \ (I -I 
-1 j ~ 2 U lO -i/ U i 



(7.50) 



Let 

M{z)={\ ~')m{z)[\ ^J. (7.51) 



Then, clearly 

(i) M{z) is analytic on C \ [aN.f^N]- 



m 



M+(^) = M_(z)(^Q ze[aN,PN]. (7.52) 

M{z) = I + {\z\-^) , \z\-^oo. (7.53) 



Thus, 



M{z) 



1 rdN ]2Skzllds 

y g 27ri ■'CKjV s — z 

i loo- ^iV~^ 

e ■i ^'^N-- 



(7.54) 



7^' 

7Ar 
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where 

with cut on [aM^P^] and the branch such that 7iv(c)o) = 1. Then 

-1 



' 7iv(^)+7jv^(^) 7Ar(^)-7iv^(^) ' 
2 (-2i) 

7iv(^)-7j^^(^) 7iv(^^)+7jv'(^) 
2i 2 



detM(z) = 1. 



At infinity we have that 



N 



At the origin, 



7^(+.0) = ( ^-^^^ = ^tan ^(1 - C) e--/^ + OiN-% (7.58) 



hence 

—p + iq p + iq 



where 



V2 
4 

We have the conjugation condition 



P,q 4 



tan^(l + C)± Jtan^(l-C) 



(7.56) 



M(+zO) = ( P P - 'A + 0(Ar-2), ^7 59) 



(7.60) 



M{z) = asM{z)a3, (7.61) 

hence 

M(-zO) = 7 + 0(iV-2). (7.62) 

\p + tq p — tq J 

8. Parametrix at the edge points 

We consider small disks D^Pn,^), D{aiy,r) of radius r > 0, centered at the edge points, 
and we look for a local parametrix P defined on D{j3N, r) U D{aN, r) such that 

(i) P{z) is analytic on {D{[3n, r) U D{aN, r)) \ (M U Sjv), where T,n = ^ ^iv 
boundary of the lenses, see Figure 6. 

(ii) P+{z)=P4z)Js{z), ze{D{(3N,r)UD{aN,r))n{RUJ:). 

(iii) P{z) = {1 + (AT-i)) M(z), 2 e dD{/3N, r) U aD(a7v, r), ^ oo. 

We consider the right edge point /Sat in detail. Note that by (|3.17|) . we have that for z G 
D{PN,r), 

-9n{z) + ^^^^^ + ^ = 1 ^ hN{fJ^)\/{fi - aN){fi - 13 n) dfi 

_ 2 
~ 3 



Jf3M (8.1) 

aN{z){z - (3n?/^, z G D{l3N,r) \ [aN,(3N], 
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where a^iz) is an analytic function in D{Pisf,r) such that 



]^ 2 



+ 0{N-') > 



Define the analytic function, 
Xn{z) ' ^ 



2/3 



so that X'j^{Pn) = o-NiPN^^^ > 0, and consider the conformal mapping, 

We will assume that the contours are chosen in D{PN,r) in such a way that 

, r 27r 
Aat : T^j^ ^ <z : argz = ± — 

Let us transform the RH problem on the matrix S{z) in D{f3i\i,r). Let 

Lemma 8.1. ^{z) satisfies the jump condition 

$+(z) = J$, 

where 

( ^^ ^\ 

for arg2; = 



J$ = < 



1 1 

1 

1 

1 1 

1 
-1 

1 

1 1 



for z G S 



Ni 



for arg z = tt, 
for z G 



(8.2) 

(8.3) 
(8.4) 
(8.5) 

(8.6) 
(8.7) 



We will use a model solution to ()8.7|1 . which is constructed explicitly in a standard way 
out of the Airy functions. The Airy function Ai(z) solves the equation y" = zy and for any 
e > 0, in the sector vr + e < arg 2; < n — e, it has the asymptotics as 2; — > 00, 

1 



2v^ 



-1/4 -223/2 



1- -^-3/2 ^^-3 + 0(^-9/2) 
48 4608 ^ ^ 



Ai'(;2) = L^l/4e-|2^/^ 1 + 1^-3/2 

^ ' 2v^ V 48 



455 
4608 



Z-' + 0(z~^/'' 



(8.9) 



The functions Ai(ci;2;), Ai(u;^2;), where uj = e 3 ^ also solve the equation y" = zy, and we 
have the linear relation, 



Write 



Ai{z) + ujAi{ujz) + uj'^Ai{uj^z) = 0. 
yo{z) = Ai(z), yi{z) = uAi{ujz), y2{z) = u'^ Ai{uj^ z) , 



(8.10) 
(8.11) 
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and we use these functions to define 
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yo{z) -y2{z) 
y'oiz) -y'2{z) 

-yi{z) -y2{z 

-y'liz) -y'^{z 

-y2{z) yi{z) 
-y'2{z) y'^{z) 

yo{z) yi{z 

y'oi^) y'ii.z 



for < axgz < 27r/3, 
, for 27r/3 < axgz < n, 

for — TT < aigz < — 27r/3, 
for — 27r/3 < aigz < 0. 



(8.12) 



Then in the sector < a.Tgz < 27r/3, 

^z) 

where 



'^'11(2;) $12(2;) 

^2l(z) $22(2;) 



(8.13) 



<fii(^) 
$12(2;) 
$21(2:) 

$22(^) 



2v^ 

27^' 
1 



_l/4 -2^3/2 

z ' e 3 



5 

48' 



2v^ 



_l/4 - 223/2 

z ' e 'i 



l/4^-|(..2^)3/2 
7 



5 

48' 



0{z 



-3^ 



48 ^ ' 



(8.14) 



2v^ 



1 + — (cj'^) 
48^ ^ 



-3/2 



+ 0(2- 



where for z z^^^, and z^^^ the principal branches are taken, with the cut on (— cxd,0). 
Since tu^ = and < argz < we have that aigcu'^z = argz — hence (cu'^zY^'^ = 

e-fzy\ {u?z)-^l^ = efz-^/\ {u?zfl^ = -z^/^, and {cu^z)-^/^ = -z-^l^. Substituting 
these expressions into (j8.13p . we obtain that 



2A 



-0-3/4 



1 i 
-1 i 



+ 



48 



-5 hi 
-7 -li 



-3/2 



+ 0(2- 



e 3^ "-^ 



Note that $(-2) satisfies the jump condition $+(2;) = $-(-2) J$. Define 
P{z) = E(;.)iV^-<|.(iV2/3A^(z))e^(-^-(^)+'^+^)- 



(8.15) 



(8.16) 



where E{z) is an analytic prefactor that has to be chosen to satisfy the matching condition 
P{z) = {I + {N~^)) M{z) on the boundary of D{(3n, r). Then 



E{z) = ^/^M{z) 



1 -1 

—i —i 



'Mz)X]i\z) 

rN\z)x-^'^\z)^ 



(8.17) 
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Recall the definition of — ( fzf^ ) ^-^^ note that 

Therefore E{z) is indeed an analytic function in D{f3N,r). 

A similar construction works for a parametrix P around the other edge point. Namely, by 
(jSUHl), we have that for z G D{aN,r), 

Yj^(z) l]\f 1 f"''^ 

-9Niz) + ^^ h Y + vrisgn(Im2;) = -y hN{ii)^J{aN - Ij){I3n - jj) dfi 

2 

= -aN{z){aN - zf^, z G /^(a^v, r) \ [oiv, /^at], 
where aj^{z^ is an analytic function in D{aj^^r) such that 

(^n{oln) = I h{aN)VPN-aN = 7 , L + 0{N-^) > 0. (8.19) 

2, [—a)y/p — a 

Define the analytic function, 



Xn{z) 



2 ( -9n{z) + + — + m sgn (Im z, 



2/3 



aN{zf'^{aN - z), (8.20) 



so that A'^(aAr) = — aiv(aiv)^/^ < 0, and then define P{z) by the formula, 

P{Z) = a3^(z)iVi'^3$(Jv2/3;^^^^))gA.(-3ivW+^ + %^)-3^^^ ^g_2l) 

where 

E{z) = V^a,M{z)a, (^^^ {XN{z)r^' 

7Ar(z)Ajv^/'^(; 



(8.22) 



-2 — Z 

Observe that the function 

j],\z)\]i\z) = - zy/'Mz))'/'. 

is analytic in D(ajv,r), hence -£^(2;) is analytic as well. 

9. Approximate solution to the RH problem 

Define 

jS{z)p-\z), ifzeD{aN,r)UD{(3M,r), 
R[z) = < _ (9.1) 

1 S{z)M [z), otherwise. 

Then, in D{aN,r) U D{PN,r) we have that 

R^z) = S4z)P~\z) = S.{z)Js{z)Js\z)Pz\z) = S.{z)Pz\z) 
= R-{z), 

on S+ U U (M \ \aM - r, + r]), 

R+{z) = S+{z)M-\z) = S^{z)Js{z)M-\z) = S_{z)M~\z)M{z)Js{z)M-\z) 
= R4z)M{z)Js{z)M-\z), 



(9.2) 



(9.3) 
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on [aN + r,(3N - r], 

R^iz) = S^{z)M-\z) = S^{z) ' M=\z) = S^{z)M:\z) 

= R-iz), 

and on dD{aN,r) U dD{[3jq,r) the jump matrix is 

J^{z) = Rz\z)R+{z) = P{z)S'\z)S{z)M-\z) = P{z)M-\z). (9.5) 
Introduce the contour E^j, which consists of the six arcs, 

= (-00, «^ - r) U U S+ U U U {P^ + r, 00), (9.6) 

where 

= dD{aN,r), = dD{(3N,r), = \ [D(a^, r) U D(/3jv, r)]. (9.7) 

see Fig. 7. The orientation of the arcs is shown on Fig. 7. 



Figure 7. The contour Er. 

For the sake of brevity we will denote 

= (-00, aN-r)U {(3n + r, 00). 

We have the following 



(9.8) 



Lemma 9.1. S{z) is a solution of the Riemann-Hilbert problem - \7.18^ if and only 

if R{z) is a solution of the following RH problem: 

(i) R{z) is analytic on C \ S^, 

(ii) R+{z) = R^{z)Jr{z), z G S/j, where 



Jr{z) 



M{z)Js{z)M-\z), on {dD{aN,r) U dD{f3N,r)) 

P{z)M~\z), on dD{aN,r)UdD{f3N,r). 

(iii) R{z) = I + 0{z-^) , z^oo. 

We evaluate the jump matrix Jr on different pieces of T,r. 



(9.9) 
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Jump matrix Jr{z) on S^. We have on dD{P]\f,r) that 

Jr{z) = P{z)M-\z) = E(;z)Ari-3$(Ar2/3A^(^))e^(-sivW+^+%^)-3]vf-i(2) 



= ^M{z) 

'-5 hi ' 
-7 -li 



1 i 
-i i 



48 



X 



X'^"\z) + O {N-^) 



= M{z) 



I + —N 
48 



6z 



^;)A-/^(.)+o(iv-) 



M-\z) 



= I + N-^j}i{z) + 0{N-^) 



where 



5{z - a) + 7{z - (3) i[5{z - a) + 7{z - 



" 96(2-/5)2(2-a)i/2a(z) V2[5(2; - a) + 7{z - (3)] 5{z -a)- 7{z - f3) 



and 



a{z) 



h{s)\/ {s — a){s — (3) ds, 



4{z - 

where h{s) is defined in fl4.16p . 

Jump matrix Jb{z) on S^. Similarly, on dD{aN it), 



(9.11) 
(9.12) 



Jr{z) = M{z) 



1 ,,-1 f 1 -Qt\ ,-3/2 



-6i -1 



A-^/^(.) + 0(iV-) 



M^\z) 



I + N-^j}i{z)+0{N-^) 



(9.13) 



where 

4(^) = 



, . 1 / 7{a -z)~ - z) i\-7{a - z) - 5(/5 - z] 

'^""^ ~ 96(a - z)\(3 - zy/^a{z) \t[-7{a - z) - 5(/? - z)\ -7(« - z) ^ 5(/? - z) 



and 



a{z) 



h{s)^/ {a — s){P — s) ds, 



4{a - z)3/2 

where h{s) is defined in (I4.17p . 

Jump matrix Jr{z) on S^. By and (ITTTIl . on S+, 

Jk(z) = / + JM, 

where 

J° (2) = e-^^'^(")M(z)aH_M(2)-\ z e E+ 

and 

^0 0' 



1 



1 





(9.14) 
(9.15) 

(9.16) 
(9.17) 
(9.18) 



From (j7.26p . (j7.31|) . we obtain that there exist constants C, 7, c > such that 

11^^(^)11 < CAT-^e""^!^'""!, (9.19) 
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where || is the sum of absolute values of the matrix elements of Jji{z). On T,]^ we also 

have equation ()9.16|) with estimate ()9.19|) . 

Jump matrix Jji{z) on S*^. By ()7.17|) . Jr{x) = I + J]i{x), where 

J°(z) = e^[^'^+W+9^-W-^^(^)-'^]M(z)a_M(2)-\ z e S^. (9.20) 

In this case, there exist C,c> such that 

||J^(^)||<Ce-^l^l (9.21) 

Solution of the RH problem for R by perturbation theory. The estimates above 
show that J'fiiz) — > as — oo. We can apply the following general result. 

Proposition 9.2. Assume thatv{X), A G solves the equation 

„(A) = / 1 /■ '^MW,,, AeE«. (9.22) 

where A- means X — iO, the value of the limit from the minus side, and Jr = I + J'^. Then 

R{z)=I-—[ 'M3Ltld^^, zeC\^n, (9.23) 

solves the following RH problem: 

(i) R{z) is analytic on C \ S/j, 

(ii) i?+(A)=i?_(A)J«(A), AgSk, 

(iii) R{z) =1 + (z-^) , z^oo. 

Proof From ^TT^ . 

R4X)=v{X), AG Sr. (9.24) 

By the jump property of the Cauchy transform, 

R4X) - R4X) = v{X)r^{X) = i?-(A) J^(A), (9.25) 

hence i?+(A) =_R_(A) Jr(A). From X^hT^ . R{z) =1 + 0{z-^). Proposition El is proved. 
Equation ()9.22|) can be solved by perturbation theory, so that 

oo 

viX) = I + J2M^), (9-26) 

k=l 

where for k > 1, 

..(A) = ~/ d,, AeE„, (9.27) 

and fo(A) = /. Series fj9.26|) is estimated from above by a convergent geometric series, so it 
is absolutely convergent. Observe that 

viW = -:^f ^^d^i, AgSr, (9.28) 

The function R{z) is given then as 

oo 

Riz)=I + J2Mz), (9-29) 

k=l 
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where 



In particular, 



1 /■ "..-MrM ^ 

Riiz) = -^ [ ^-^d^j.. (9.31) 
27rz Jt.„ z- n 



10. Large N asymptotics of the recurrent coefficients 
From ()6.9p . ()7.H) and ()7.15p . we obtain the formulae for the recurrent coefficients: 

hNN = -27ne^'HSi)i2, 

RnN = ('S'l)2l(5'i)i2, 
QnN = TTTT^ + ('S'l)ll. 



where 



where 
Therefore, 

By (unzD, 

hence 



^ ^ j°(A)rfA--i. / v^{X)J%{X)dX-... 



2wi y ' 2m 



We will call the first term on the right. 



a linear term, the second one. 



(10.1) 



S{z) = I+^ + ^ + 0{z-'), z^oo. (10.2) 

z z^ 



By dHH), S{z) = R{z)M{z) for large z, hence 

5^ = Mi + i?i, (10.3) 



M(z) =/+ — + 0(;2-2), R(^z) = I + — + 0{z-^). (10.4) 
i?^^ = (Ml + i?i)2i(Mi + i?i)i2. (10.5) 



i?^^ = (^^^^^ + ^^L^[{R,)u - {RM + (i?l)l2(i?l)21. (10.7) 

By 

R, = -± [ viX)J%iX)dX 



(10.^ 



(10.9) 



R? = I v,{X)J],{X)dX, (10.10) 
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a quadratic term, etc. By definition, we have that 

i?i = i?« + i?f) + ... (10.11) 

First we evaluate the hnear term. 

Evaluation of the linear term. Denote 

R1 = -^l JrWcIX, a = a, /?,+,-, oo, (10.12) 

so that 

i^Jl) =R'^ + Rf + R+ + R-+ R^, (10.13) 

Let us evaluate R^, Rf, R^, and R'^. 
Evaluation of R^, Rf. By ^TW^ . 

Rf = f 4W dz + 0(iV-2), (10.14) 

which can be evaluated by taking the residue at z = p. The result is 

Rf = -N^'— f .,ff + " , - + 0{N-'). (10.15) 

192 \*(ll/3 - a) -3/5 -a J 

A similar expression holds for R^. Namely, 

1 / -3a -p i{lla-p) 



R^ = -N-' , ] +0(N-'). (10.16) 

By taking into account terms of the order of A^~^ in ()9.10j) . we obtain the error terms in 
pO.lSp . p0.16p as N~'^Ca,i3 + 0{N~^), where Ca^p are some explicit matrices. 

Evaluation of Rf. In the usual case of a random matrix model with an analytic potential 
V{M) independent of A^, the terms Rf, which represent the integral over the lenses boundary, 
are exponentially small as ^ oo, see It is not the case in our situation because of the 
series of poles of the function e"^*^^*^^) on the imaginary axis. By (j9.17|) . 

Rt = — - I JU^)dX = — - [ e-^^'^WM(A)a+M(A)-MA. (10.17) 
27Ti Jy,+ 27ii Jy,+ 

From ()7.27|1 we obtain that the function e"^*^^*^^-* has simple poles at the points 

tN-'y„ y^. = -^, J = 1,2,..., (10.18) 

27 



Z = Zj 



and by (j7.34|l . the residue at Zj is equal to 



ti — 1 ) sin V ' ~~~ 

Res [e-^^^(^)] = e'^'^kNiVj)—^ ^(1 + O(iV-i)), (10.19) 



iV ( JL _ 1 



By using ()7.36p we reduce this to 

Res [e-^^^(")] = iCje'^'^N-^^il + 0{N~')), (10.20) 
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where 

= 1 + ^, C, = ^e-fe)(-iysin (-^) . (10.21) 

Observe that > 1. From ()7.46p we obtain that 

Res le^^'^^'^ = -iCje~'^'^N-''^{l + 0{N-^)), (10.22) 

Z = -Zj 

Let us deform the contour up, crossing the poles. Every time we cross a pole, the 
residue at the pole appears on the right of (jlO.lTj) . while the integral becomes smaller than 
the contribution from the pole. This gives the asymptotic expansion as ^ oo, 

oo 

~ - y Res re-^^^(^)] M{z,)a+M{zj)-\ (10.23) 

' ^ Z = Z'i 

i=i 

where the j-th term is of the order of N~'^k For our purposes it will be sufficient to consider 
terms with /tj < 2 only, 

= - y Res [e-^^'^(^)] M{zj)a+M{zj)-^ + 0{N~^-'), (10.24) 

where 



Z = Zj 

j: Kj<2 
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Jo 

In fact, since Zj = 0{N~^), we can replace M{zj) by M(+iO), 



10.25) 



i?+ = - y Res [e-^^^^")] M(+zO)a+M(+iO)^i + 0(A^-2-"). (10.26) 

^ Z = Zj 

j: Kj<2 



(10.27) 



Let us rewrite this in terms of the matrix elements, 

(i?+)i2 = M+'4 + 0(iV-2-^), (RtU = -M+\4, + 0(iV-2-), 
{Rt)n = -M+M22 Js^i + 0{N-^~'), (i?+)22 = M+M224 + 0(iV-2-^), 

where 

J+= J2 Res [e"^^'^^")] = ze^^^ ^ CjN"''^ + 0{N-^~'), (10.28) 

j: K,j<2 ^ j: K,j<2 

and are the matrix elements of the matrix M{+iO). By applying ()7.59p . we obtain that 

iRt)i2 = ip- iqfJti + 0(iV-2-), (i?+)2i = -{P + ^qfJti + 0(iV-2-), (10.29) 
Similarly, we evaluate the contributions from the contour as 

(i?r)i2 = (P + iq?J2i + 0(iV-2-^), (/?r)2i = -(P - ^g) V2-1 + 0(iV-2--), (10.30) 
where 

J2-i = - y Res [e^^'v^")] = ze"^^" ^ C^-A^-''^- + 0(A^-2-£)_ ^^0.31) 

j: Kj<2 Kj<2 
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By combining ()l().29j) and (jlO.HOj) . we obtain that 

(i?+)i2 + (i?r)i2 = «2 [{p^ - q^) cosiNu) + 2pqsm{Nu)] ^ C^N'^^ + 0{N- 

j: Kj<2 

{Rt)2i + (i?r)2i = -^2 [{p" - q^) cos{Nu) - 2pqsm{Nuj)] ^ CjN~''^ + 0{N- 

j: Kj<2 
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2-e\ 



2~e\ 



From ()7.60|) we find that 



2 2 1 o 1 7'"C 

P - g = 2 ' 2pg = - tan y 



(10.32) 
(10.33) 



hence 



(i?+)i2 + (i?r)i2 = ^ 
(i?+)2i + (i?r)2i = 



cos(A^c<j) + tan ( ] sin(A^u;) 



■2-e\ 



j: Kj<2 



cos{Nuj) — tan ( -^^ J sin(iVu;) 



J2 CjN^'^^ + 0{N- 

j: Kj<2 



2-e\ 



(10.34) 

Evaluation of i?^. From (j9.2H) we obtain that R'^ is exponentially small as — oo, 

= 0{e-'">^). (10.35) 

Summary for the linear term. The evaluation of the linear term can be summarized 
as follows: 



{R\ 



(1). 



112 



-N- 



5i{(3 — a) 
96 



+ i 



(i?f))2i = -N 



96 



— t 



cos{Nu!) + tan ( ] sin(A^co' 



cos(A^co') — tan ( — ) sin(A^Ci;) 



j: Kj<2 



J2 C,N~^^ 

j: Kj<2 



(10.36) 



-2-e\ 



where ^ c"2i some constants. 



Evaluation of the quadratic term. We obtain from ()9.28|) . ()10.10|) . that the quadratic 
term is equal to 



R\ 



(2) 



dfidX. 



We can split it as 



where 



i?f)= J2 ^i''' A = {a,P,+,-,oo}, 



a,beA 



R1'' 



i2ny 



A_ — /i 



dfidX. 



;i0.37) 
;i0.38) 

;i0.39) 



40 



PAVEL M. BLEHER AND VLADIMIR V. FOKIN 



If a ^ b then we can replace A_ by A and in this case we obtain that 



(27r)^ 



dfidX, a ^ h. 



(10.40) 



It is tempting to say that R{'^ = —R^' , but in general it is not true, because the matrices 
J^(A) and J^(/i) do not commute. By ^TK^ J'^{z) is analytic on and by (FT^ it is 
exponentially small in A^l^;], hence i?"''' is exponentially small in N, if at least one of a, b is 
equal to oo, 

iRfl < Coe-^«^ 
From (|9.10p we obtain that 
1 



Co, Co > 0; a = oo or b = oo. 



2t\% JyP A- — [i A'=/3 A — p 



hence 



Similarly, 



^ = Ar-2Res Res 



_ ;Y-2p^g Res 



li=a 



4(A) 
A — a 



-3\ 



+ 0(iV 



+ 0(A^-^) 



The cross terms are evaluated as 



(10.41) 



(10.42) 



;i0.43) 



10.44) 



'- Res [J'M] Res [4(A)] + 0(iV-3) 

P — a M=/3 A=a 

'-^Res [4(/i)] Res [4(A)] + 0(iV-=') . 
p — a A=/3 



Thus, 
where 



ci = Res [4(A)] Res 



X — a 



+ Res [4(A)] Res 



4(A) 
A-/? 



+ ^T^Res [4(A)] Res [4(A)] - -r^Res [4(A)] Res [4(A)] 



Let us evaluate R^'^ . Consider 



vtiX) 



27ri Jy,+ A_ — yU 



;i0.45) 



(10.46) 



(10.47) 



'10.48) 



By deforming the contour of integration up, we obtain the asymptotic expansion of fj^(A) 
as iV — cxD, 



oo ^ 

+ (A) ~ -$^3^-^Res [e-^^'^(^)] M{zj)a+M{zj)-\ 



(10.49) 



'3 -J 
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Now we substitute this asymptotic expansion into the formula, 

RV^ = -7^. [ vt{x)r^{x)dx, (10.50) 

and move the contour of integration up. This gives the asymptotic series, 

J2 Res [e-""^^^'^] Res [e-^'^^'^'^M{z,)a+M{zj)-^M{zk) 

Zj^ Z n Zj Z Zj^ 

(10.51) 



'-^-NGn(z) 

X a+M[Zkr^ + 2^Kes ye~"'^"'-^'\ Kes 



Z-Zj 



M{z,)alM{z,)-\ 



Observe that the last sum is equal to 0, because cr^ = 0. Furthermore, since M{zj) = 
M{+iO) + 0{N~-^), we obtain that 

M{zj)-^M{zk) =1 + 0{N-^). (10.52) 

When we substitute / for M{zj)~^M{zk) in the first sum in p0.5ip . we again get 0. When 
we substitute 0{N~^) for M{zj)~^M{zk), we get a term of the order of 0{N^'^'^'^). Thus, 

R+'+ = 0(A^"2"^). (10.53) 

Similarly, 

R-'- = 0(A^"2"i). (10.54) 

Observe that by (jlO.'^ljl . 

= 1 + > 1- (10-55) 

TT — Z7 

Consider now 



Ri=-7r- / vtiX)JnWdX. (10.56) 

When we substitute asymptotic expansion ()10.49|) into this formula and move the contour 
of integration, S^, down, crossing the poles of J]j(A), we obtain the asymptotic expansion, 

~ y ^^Res re-^^-(^)l Res [e^^-^^)] 

j^i^i + ^fc"="^ ''=-'^' (10.57) 

X M{z,)a+M{zjy'M{-Zk)a+M{-Zk)-\ 

Since 

M{zj) = M(+zO) + O(Ar-i), M{-Zk) = M(-zO) + 0{N-^), (10.58) 
we have that 

Miz,)a+Miz,)-'Mi-Zk)a+Mi~Zk)-' 
= M{+iO)(T+M{+iO)-^M{-iO)(T+M{-iO)~^ +0{N~^). 

From (|7.56p we obtain, by a direct computation, that 

a+M{+iOy^M{-iO)a+ = -a+, (10.60) 
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and from ^SM, ^SM, that 

Res [e-^^^(")] Res [e^^^^''^] = CjCkN-''^~^''{l + OiN-^)). (10.61) 

Hence 



oo _ 

R-'^ N-^i-^^ [CjCkM{+iO)a+M{-iO)-^ + 0{N-^)] . (10.62) 

j,k=i + ^'^ 

Observe that 

a+M{-iO)-^M{+iO)a+ = a+, (10.63) 
and, therefore, a similar computation for Ri~ gives that 



~ — ^ — iV-'^^-«'= [CjCkM{-iO)a+M{+iO)-^ + 0{N-^)] . (10.64) 
j,k=i + 

Since 

M{-iO)a+M{+iO)~^ - M{+iO)cr+M{~iO)-^ = I, (10.65) 

we obtain that 

j,fc=i + ^'^ 

If we restrict this matrix formula to the elements (12) and (21), then we obtain that 

(i?+'")l2 + (i?r+)l2 = 0(iV-2-i), (i?+'-)2i + (i?-'+)2i = 0(iV-2«i), (10.67) 

because 7-^^ = 0{N). Finally, the cross terms of the form i?^'^, where a = ±, 6 = a, /?, or 
vice versa, are estimated as 

Rf = 0{N-^-^'), a = ±, b = a,P, or a = a, /3, 6 = ±. (10.68) 

Summary for the quadratic term. By combining formulae ()10.4H) . ()10.46|) . p0.53|) . 
(|l().54l) . (Il().67|) . and (ll().68j) . we obtain that 

{R?)i2 = (ci)i2iV-^ + 0(iV-2-), (i?f ))2i = (ci)2iiV-2 + 0(iV-2-), (10.69) 

where the matrix ci is given in (jl0.47|) and e > 0. 

Evaluation of the higher order terms. The higher order terms, R'"i \ k > 3, are 
evaluated in the same way as the quadratic terms, and we obtain that 

{r[%2, {R{%i = 0{N-'-n, k > 3. (10.70) 
Consider, for instance, the cubic term, 

^^^^Jsn V 27r?y Js^ Je^ Js^ - /^)(/^- - '^J 

(10.71) 

As for the quadratic term, we split into a sum of terms i?"'*'^, and the only nontrivial 
terms in regard to estimate p().7()|l are Ri'~'~^ and R^'~^~. We have that 
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We move the contour of integration up and the one down, and obtain the asymptotic 
series, 

~ y J .Res re-^^-(^)l Res [e^^-^^)] Res [e-^^-^^)] 

jfji^ i^j + ^k){Zk + Zl) z=z, z=-z^ z=zi (10.73) 

X M{zj)a+M{zj)-^M{-Zk)(y+M{-Zky^M{zi)a+M{ziy^. 
By using (jl().58|) and (jlO.eiOjl . we obtain that 

oo 

~ I V ^ N--^--^--^ [C,CfcQM(+^0)aH_M(-zO)-i + 0{N-^)\ . 

^k){Zk + Zl) 

(10.74) 

A similar computation for K^'^'^ gives that 

CO ^ 

^r^'"~-^ E ^ ^7 ^iV-'^^-''''-''' [C,CfcQM(-z0)a+M(+^0)-i + O(iV-i)], 

(10.75) 

and by using (|10.65p . we obtain that 



oo ^ 

hence 

(i?|'-'+)i2 + (i?r^'li2 = ©(AT-^'^^+i), (i?+'-'+)2i + (i?r'^'l2i = 0(Ar-3-+i). 

Since 3ki — 1 > 2, we obtain estimate p0.70j) for R^'''^ + i?^'"*"'". It is straightforward to 
get the estimate, 

^a,6,c ^ (9(^^-2-.)^ (lQ_7g) 

for all other combinations of a, 6, c and hence (jl0.70p follows. The same argument holds for 
> 3. 

Evaluation of Rnn- Let us go back now to formula (jl0.7|) and evaluate the terms on 
the right in this formula with an error term of the order of N~'^~^. From ()5.10p . 

PN-aN = + N-^- ^^^^ + 0(Ar^=^), (10.79) 



cos 



2 

hence 



^ 3(7r — 27)cos^ 



2cosfy 6(7r - 27) cos2 ^ 

Next, from (ITIOHll . (fTTTHHl) . and (ll().7()jl we obtain that 

^~"[(i?i)i2-(i?i)2i] = ^-^cos(iVa;) 5^ C,N-^^ +cN-^ + OiN-^~' 



Ai 

j: Kj<2 



cos(iVcj) CjN-"^ + c°A^-2 + 0{N' 

j: Kj<2 



■2-e\ 



(10.81) 
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where 

c, = = ^^i-iy sin , (10.82) 

^ cos — ^~ — 

and c° is a computable constant. From ()10.7|1 we obtain now that 

Rnn = [ — ^ ) + cos(iVcu) J2 Cj-^""' + c^"^ + ^(A^"^"') , (10.83) 

y2C0S^y j;Kj<2 



where 



+ Co. (10.84) 



6(7r - 27) cos2 ^ 

Here the first term in the expression for c comes from the difference ( ^^~"^ )^ — (^^)^, see 
(110. 80|) . while the second term, cq, is determined by calculations of other terms of the order 
of iV~2 on the right in formula ()10.7|) . The constant cq can be evaluated explicitly by tracing 
down all the terms of the order of A'^"^ in the above computations. To avoid these somewhat 
tedious computations, we will use the fact that we know the exact expression for Rnn on 
the free fermion line. 

Observe that Cq is calculated in terms of contour integrals around the turning points un 
and Pn, and it depends only on the limiting values of the end points, a, j3. The exact values 
of a, /3 are given in (j4.5|) and they depend on the parameter C, only. This implies that cq is a 
function of the parameter ^ as well, Cq = Co(C), and it is independent of 7. To find an exact 
value of co(C), consider the free fermion line 7 = f . In this case c = 0, which gives 

Co(C) = -— (10.85) 
48 cos^ -Y 



Thus, 



717^ TT^ 



(10.86) 

6(7r - 27) cos2 f 48cos2f 



2 

This proves formula ()2.10|) and hence Theorem 11.11 

11. Proof of Theorems 11.31 and 11.41 

We omit the proof of Theorem II. 2^ because it follows from Theorem 11.41 
Proof of Theorem 01 By (fT3H|) and (IT^ . 



Q'^Fn Rn 1 



R + cos{Nuj) CjN-'^^ + cN-'^ + 0{N' 

j: Kj<2 



It is easy to check that 

d^F _ R 

hence ()1.57j) follows. Theorem 11.31 is proved. 
Proof of Theorem 01 By 

n^R , 



;ii.i) 



(11.2) 



Rn = ^e"", (11.3) 
7 
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where 



bn = cos{nuj) djU '^^ — nn ^ + 0( 

j: Kj<2 

and = ^ , K = — -|. From ()1.27j) . ()l.H2j) we obtain that 

c 1 



n 



K = —- 



R 12 37r(7r - 27) 
From ()1.25|) and pi.3|) we obtain that 

N(N-l) /JV_1 ^ 2 



' (n 



where 

Bn = {N- l)bi + {N- 2)62 + ■ ■ ■ + bN-i, 

hence by fOlj) . 



= iV-^ln- ^^^^^ = iV- 



, N{N-1)^ R 

N\nho + ^- ^In — 

2 7"^ 



-In— + CoAr-' + Ar-'5^, 



2 72 

where Co is a constant. Let us evaluate B^. We have that 

Bn = N{bi + 62 + ■ ■ ■ + &iv) - &i - 262 NbN 

and 

00 

61 + 62 + ■ ■ ■ + = 5 - ^ 

n=JV+l 

where 

CXD 

n=l 

It follows from (fTOj) . that 

00 

J2 b^ = -KN-' + 0{N-'-n, 

n=N+l 

because 

00 

n-"^ cos{nu) = 0{N-^^), 0<u}<2n. 

n=N+l 

It also follows from ()11.4|) . that 

N 
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where Ci is a constant, because 

N 

J2 ^""'^^ cos(ncj) = C{kj) + 0{N-^^+^), 0<u<27i. (11.15) 

n=l 

Thus, 

Bn = C2N + KlnN + C3 + 0{N~'), (11.16) 
where C2, C3 are some constants, hence from ()11.8|) we obtain that 

Fj^ = F + cqN-^ + kN-"^ \nN + C^N-'^ + O^N-"^-'), (11.17) 
where Cq is a constant. This imphes that 

Zjv = Ce^'/+^^o;vK (1 + 0(Ar^^)) , (11.18) 
where C = e*"^. To finish the proof of Theorem II. 4[ it remains to prove the following lemma. 
Lemma 11.1. cq = 0. 
Proof. By 

\ 2n+l 



K = Kn, (11.19) 

and by ^UU} . 

/i„„ = -27r«e"'"(5i)i2. (11.20) 

Observe that by ()5.(i7|l . 

/^ = / + 0(^-2), / = 21n(/3-a) -2-41n2, (11-21) 

and by (fTn3|l . (ITin)|l . 

(5i)i2 = (Mi)i2 + (i?i)i2 = (1 + 0{n'')) . (11.22) 

Therefore, 

K = + o(n-i)) , (11.23) 

and 

A'-l /N-1 \ /in ^\ N-1 



tn = 

n=0 \n=l 



X exp (^ ^^^^ + O(lniV) ) . 



11.24) 



By applying (111.211) . (OH|) . and (gSl), we obtain that 

tn = CNexp{N^F + 0{\nN)) , (11.25) 
where Cat does not depend on 7 and t. By fll.7|l and ()1.37j) . this implies that 

Zn = CNexp{N^f + 0{\nN)) , (11.26) 

where Cn also does not depend on 7 and t. Since on the free fermion line, = 1 and 
/ = 0, we obtain that 

InC^v = O(lnAr), (11.27) 
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hence 

Zn = exp{N^f + 0{\nN)) , (11.28) 
so that Co = 0. Lemma fl 1.11 is proved. 

Appendix A. Large asymptotics of A{N) and A{N; 3) 
Large N asymptotics of A{N). We will find in this appendix the large asymptotics 

of 

^ ^ ii (2n)!(2n+l)! ^ ^ 

We start with the asymptotics of 

7V-1 

a{N) = Yl n\. (A.2) 

n=l 

We have that 

N N N 

In a{N) = ^(A^ - n) In n = ^ (^1 - ^) (in ^) A^"^ + ^(A^ - n) In A^. (A.3) 

n=l n=l n=l 

In addition, 

N 



ri=l 

where ({s) is the Riemann zeta-function. This gives 

, ,,,, NHnN 3N^ N\n{27T) InA^ 1 ^ ^ 

In a(A^) = + ^ — + C (-1) n + (A.5) 

V ; 2 4 2 12 ^ ^ ^ 240A^2 ' v ; 



so that 

Consider now 

We have that 
where 

Observe that 



N-l 



a{N) = Y[n\ = N^e-'i^\2Ti)^N-T-2e^'^-^^-^+-. (A.6) 



n=l 



N-l 



n=l 

lna3i(Ar) = hi{N) + 6o(A^) + &-i(A^) (A.8) 

N 

b,{N) = - n) ln(3n + j), J = 1, 0, -1. (A.9) 



n=l 



6o(Ar) = ^(AT - n) \n{3n) = (^^^)^(^ — 11 + \na{N), (A.IO) 



n=l 
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hence by (IA.5|1 . 

(ln3)iV(iV-l) , NHnN 3N' , iVln(27r) IniV 1 , 
2 + ^ r + ^ 12" + ^ 240A^ + - • • 

Now, 



W{N) + 6_i(iV) - 26o(iV) = f^(iV - n) In (l - ^] 

n=l ^ ^ 



, 3V3\ InAT 2 

A^ln — ^ H h 7o H ^ + 

' 27r / 9 ^ 243Ar2 ^ 



where 70 is a constant, 

TV 



7o = lim 

Af^oo 



— n In I 1 

n=l ^ 



1 \ In AT 
9^/ 9~ 



Therefore, 



,,,, 3(ln3)A^(A^- 1) 3A^2in^ 9;y2 3Arin(27r) InA^ o./. .^ 

1 /^3V3\ InAT 2 

+ • • • + A^ln -— + — — + 7o + 7— 7T7TT + • • • , 



and 
Finally, 



80A^2 V / ^ ^^^^^ 

TV-l 2Ar-l 

JJ [(2n)!(2n + 1)!] = JJ n! = a{2N). 



n=0 n=0 



By (jSJ 



a(2iV) = (2iV)2^'e-=^^'(27r)^(2iV)-i^e^'(-^)-^+-. 
Thus, flA.l|l reduces to 

^^^^ _ a3i(A^)a(A^) _ A^— 3 — 6-4^ (27r)TA^-36 A^— (27r)TA^- 



a{2N) {2Ny^'e-^^\27c)^{2N)-T2 



g3C'(-l)+7o-^5j§^+...gC'{-l)-^+. 

X 



C'(—l)— + 



By simplifying, we obtain that 



Ar2 



3^3 \ A 115 , ^,,,_3 



where 

C = 2^e3«'(-^)+^°. 
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Large asymptotics of A(N; 3). From ()1.48j) we have that 



A(2m;3) = 3' 



,2 m! 



m— 1 

n 



(3A; + 2)! 



m— 1 



A(2m + 1;3) = 3™'+™ JJ 



fc=0 



(3fc + 2)! 
{m + k + 1)! 



n 2 



cf. [?]. Let us start with A (2m; 3). We can rewrite it as 



A{2m; 3) = 3' 



,2 m! 



where 



Observe that 



(3m)! 



m—l 



a32{m)a[m) 
a(2m) 



«32(m) = n(3^ + 2) 



fc=0 



m— 1 



a32{m) = a3i{m) J]^ (3A; + 2) = a3i(m)3' 



k=0 



r (.n + 1) 
r(|) 



hence from ()A.22|1 and ()A.18|) we obtain that 



A{2m; 3) = 3' 



,2 m! 



(3m)! 



3"T (m + I) A{m) 



We have that 



and 



(3m)! V3 

r(m + i) 



■ g 9m 



1 _ 1 _L 1 



777, 3 e 9™ 162m2 



m! 



By combining this with asymptotics ()A.19j) . we obtain that 



A{2m; 3) = C3 



where 



3 
2 



4m^ 



3-'"(2m)T8 1 + 



77 



7776m2 



Oim 



2fn 



,6C'(-l)+27o 



Consider now A(2m + 1; 3). From (jX^H), 



N ™ (3m)!m! , , 
A{2m + 1; 3) = 3" ^..^ \ A{2m; 3). 



[(2m)!] 



By using the Stirhng formula we obtain that 



(3m) !m! 
[(2m)!]2 



4m 
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2m 



g 36m 144m^ 



2m + 1 

By combining these formulae with ()A.28|) . we get 

(V3)-(2™+i)(2m + l)T8 1 + ^ + 0{m-'^] 



A(2m+l;3) = C3(0 



We have: 



7776m^ 

Appendix B. Proof of formula (|4.1Uj) 
uj{z) dz = zuj[z) — I zuj'{z)dz. 



From (jOl . 



^ ^ = — 



1 



1 

m 



y/l3{z-a)-t^-a{z- 13) 2z 



/z~a \J z—fi 



{P — a)z 
1 ^/3(-a) / /z^' 



-a 



hence 



TT — a)z \ \ z — a y z — (3 
^/K- 



uj{z) dz = zuj{z) + 



-a 



'Kz^ {z — a){z — (3) 
dz 



vr J ^^z-a){z-f3) 



( , 

zuj[z) H — log I V -2 — a + y z — (3 



TT 



From (gSl), 

—a) = TT, hence 

^(2;) = + 21og {^^z-a + - + C. 

As 2; — > 00, 

(7(2) = logz + 0{z~^) = z\z-^ + 0(z"2)] + 2[log(2v^) + 0(z-^)] + C, 
hence C = -1 - 2 In 2, and (p3T1|l follows. 

Appendix C. Proof of Proposition 15.11 
From (jEHI)> we have that 

1 



2vr - aiv)(/3Ar - x) 



and 



27r 



QJV 



dx = 1. 
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In (Idlj) . ()(I2j) we can rewrite the integrals as the contour integrals, 



Fn(oin, Pn) — - — : 



1 



G^ioiN, Pn) — - — : 
Am 



- ajv)(^ ■ 



zdz. 



N 



(C.3) 



zdz. 



IT, \/{z- aN){z - (3n) 

where the function ^y {z — a]\[){z — (3^) is considered on the principal sheet, with a cut on 
[aNyPw], and is a positively oriented contour on the complex plane around [0^,^^], 
which consists of the two circles, {\z — a\ = e} and {\z — f3\ = e}, and the two intervals, 
[a + e, P — e], along the lower shore of the cut, and [P — e,a + e], along the upper shore, see 
Fig. 8. It follows from representation ()C.3|) that both and are analytic functions of 

«Af, Pn- 





Figure 8. The contour F 



By (ini, 

VUz) = V'{z) + f{Nz), (C.4) 

where 

V{z) = z sgnRe z — (z, V (z) = sgnRe z — (, (C.5) 

and 

TT 7r/7r \ f ^ \ 

f(z) = — cothz 1 coth^; 1 — sgnRez. (C.6) 

^^27 27 V27 ; V27 ; ^ ^ 

Therefore, we can rewrite equations ()C.1|) . ()C.2|) as 



( If Viz) 1 /■ 

F{aN,PN) = -r- (P dz = -— (b -y 

47r^ Jr. y/iz-a^)iz-Pj,) 47rz Jr. y^iz - a^)iz - P^) 

1 / zV'(z) , 1 / zf(Nz) 

G{aN,PN) = — (P dz = l-—i ■'^ > ^= 

/r, ^{z-aN){z~ Pn) 47r^ J^^ ^^z-aN){z- Pn 



(C.7) 

We will assume that aN — a = 0{N ^) and Pn — P = 0{N ^) as iV — > 00, where a and P 
solve the system 

F{a,P) = 0, 

G(a,/3) = 1, ^^-^^ 
and we will prove the existence of uniPn by using the implicit function theorem. Observe 
that a and P are given by formulae ()4.5p . 
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The function f{z) is exponentially decaying as \R.ez\ oo, and this allows us to evaluate 
the integrals on the right in ()C.7|) asymptotically, as — oo. Namely, 



1 



f{Nz) 



^J{z - aN){,z - Pn) 
1 



dz 



1 



2tiN 
fix 



NI3m 



fix) 



-NaM A/(iV - aN)if3N - N ^x) 



dx 



Observe that that /(— x) = —fix), hence 



xiajy + Pn) 
2N{-aN)pN 



(C.9) 



dx + OiN-^). 



fix)dx = 



(C.IO) 



and 



xfix)dx = 2 / X 



vr 1 



2|2 



7 e""^ - 1 
7 



TT 



7 



^(--2) 1 

e T — 1 



vr vr — 27 



u 



-du 



27^ 



(C.ll) 



3(7r -27)' 



Also we can replace oat, (3jq for a, /5 in ()C9|) and use formulae ()4.6|) . This gives us that 
1 f f{Nz) 



/r, v/(z-a^)(z-/?Ar) 
Similarly we obtain that 

zfiNz) 



r. 7^tan^ , 
dz = -N-^—^, §— + OiN-^) 



37r2(7r - 27) 



dz = -N~ 



7 



/r, v^(z-ajv)(^-/3jv) 
Thus, system fjC.7|) reduces to the following one: 



37r(7r — 27) 



+ OiN~ 



(C.12) 



(C.13) 



N) 



r. 7^tan^ , 
N-^--^, §— + 0(A^-3), 



37r2(7r - 27) 



G(ajv,/5jv) = l + iV" 



7^ 



(C.14) 



37r(7r — 27) 

In the linear approximation the latter system reads 



+ OiN- 



(ttAT - /3) + i(3N - f3)Fp^{a, (3) = N- \^'j^'' [ ^ + OiN~^' 



(ttTv - Oi)Gaj^ («,/?) + iPN - P)Gf3N io^,P) = N 



37r2(7r - 27) 

„2 



-2 



7 



(C.15) 



37r(7r — 27) 

The coefficients of this linear system can be evaluated explicitly. Namely, we have that 



F{aN, Pn) = - o + - arcsm , 

^ X Cif^N + ftAf) , \/Av(-«Jv) , /^at + Oat . /^at + Oat 

G{aN,PN) = 1 1 :: arcsm 



(C.16) 



TT 



2n 



Pn — OiN 
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which gives that 



Faj, {a,P) = ^ (^1 + sin ^ ^ , F^^ = -L _ gin 

By solving system ()C.15|) . we obtain that 

,72 (2 sin ^-1) ^ 
3{n — 27) cos -f- 

^7^ (2sin=f + 1) , 
3(7r — 27) cos -f- 



(C.17) 



(C.18) 



The determinant of system ()C.15|) is not equal to zero, and this guarantees, by the implicit 
function theorem, that there exists a solution to ()C.1|) . ()C.2|) . which has the same asymptotics 
()C.18|) . Proposition 15. II is proved. 



Appendix D. Proof of Proposition 15.21 



To prove ()5.26|) . we would like to replace rAr(yu) and r^^x) in ()5.16|) by rAr(O) and to 
estimate the error term as 0{N~'^). Fix any < r < imin{— q;,/3} . 
Case 1, G [a + r, [3 — r]. From we have that 



27r2 



27r2 



:P.V. 



' rjv(/i) 



f{Nx)dx 



X 



+ 



27r2 



f{Nx)dx 



:D.r 



]Ri\[ajv,/3iv] ^^ 



X 



Due to estimate ()5.22|) . the second integral is exponentially small as — > cxd, hence we can 
drop it. In the first integral we can drop the sign of the principal value, because the function 
under the integral is smooth, and we can restrict the limits of integration to (a + |) and 
(/3 — |) plus an exponentially small term. Finally, the function 



r7v(x) 



1 



/i — X 



is a uniformly bounded analytic function in a fixed complex neighborhood of (a;, /i) G [a + 
|, /3 — |] X [a + r, /5 — r], hence 



n 


'Up) ; 







f{Nx)dx 



/i 



X 



0{N- 



(D.2) 



because / is an odd exponentially decaying function. This proves Proposition 15.21 for fi G 

[a + r, P — r]. 
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Case 2, jj, E [oat, Pn] \ [a + r, (3 — r]. Suppose /x G [/5 — r, Pn]- From 

w X _ VMaO r^^' f{Nx)dx VMaO T""" f{Nx)dx 



27r2 iaj^ {fi - x)^yrN{x) 27r2 7^_^2r {fi - x) ^/r^^ 



(D.3) 



The first term is exponentially small as — oo (because / is exponentially decaying), and 
the second one is 0{N^^) (because / is odd and the integration is with respect to a smooth 
kernel). Let us consider the third term. We can rewrite it as 



27r2 ■ •y^_2,. {i2~x)^yrN{x) 27r2 ' ' 7^.2^ (/i - x) ^/M^O ' 

The second term is evaluated explicitly as const. /(A^/J^v), and it is exponentially small as 
^ oo. We can represent the first term as a half-sum of contour integrals over two 
contours, T±, where r+ (r_) goes from /5 — 2r to /i — 5, where 6 = ^{Pn — lAi then along 
the upper (respectively, lower) semicircle of radius 5 centered at /x, and then from /i + 5 to 
(3n- The both integrals are exponentially small as — > cxd, hence the third term in ()D.3|1 is 
exponentially small, and pjv(/^) = 0{N~^) when fi G [P — r,^^]. From ()5.24j) we obtain that 
k{Nfi) = 0{N-'^) when /i G [/3 - t/Pn]- This proves (j^:^ for /i G [/3 - r,(3N]- Similarly, it 
holds for yU G [a^v, a + r]. Proposition 15.21 is proved. 
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